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We calculate generalized parton distribution functions in a field theoretic formalism using a covari¬ 
ant Bethe-Salpeter approach for the determination of the bound-state wave function. We describe 
the procedure in an exact calculation in scalar Electrodynamics proving that the relevant correc¬ 
tions outside our scheme vanish. We extend the formalism to the Nambu-Jona-Lasinio model, a 
realistic theory of the pion. We go in both cases beyond all previous calculations and discover that 
all important features required by general physical considerations, like symmetry properties, sum 
rules and the polynomiality condition, are explicitly verified. We perform a numerical study of their 
behavior in the weak and strong coupling limits. 

PACS numbers: 24.10.Jv, 11.10.St, 13.40.Gp, 13.60.Fz 


I. INTRODUCTION 

Hard reactions provide important information for unveiling the structure of hadrons. The large virtuality, Q^, involved 
in these processes allows the factorization of the hard (perturbative) and soft (non-perturbative) contributions in their 
amplitudes. Therefore these reactions are receiving great attention by the hadronic physics community. Among the 
hard processes, the Deeply Virtual Compton Scattering (DVCS) merits to be singled out, because it can be expressed, 
in the asymptotic regime, in terms of the so called generalized parton distributions (GPDs) Q,iii The GPDs 
describe non-forward matrix elements of IDht-cone operators and therefore measure the response of the internal 
structure of the hadrons to the probes illlli There is much effort under way related to the measurement of 
these functions. 

Due to the impossibility at present to determine the GPDs from Quantum Chrom ody namics dire ctly , models have 
been used to provide estimates which should serve to guide future experiments (l3, [IB HB HB Ull • The aim of 
our work is to perform such a calculation in a field theoretic scheme which treats the bound-state in a fully covariant 
manner following the Bethe-Salpeter approach. In this way we would like to preserve all invariances of the problem. 
For simplicity we shall use mesons as initial and final states. 

We define a scheme, to calculate the electro-magnetic interaction of hadrons, which separates the soft parts, where 
we use a non perturbative treatment, and the hard parts, where the conventional perturbative treatment is applied. 
The scheme preserves gauge invariance to leading order and leading twist. 

In order to describe the soft part we start by using as a test of our ideas a model based on the (f>‘^ field theory. This 
theory has the advantage of being renormalizable and therefore any contribution can be analyzed properly. We show, 
for example, that a correction to the hand bag contribution to first order in the strong coupling constant, which goes 
beyond our scheme, vanishes. 

We next proceed to the main development of this paper, namely to perform the study of the GPDs of the pion by 
using the Nambu-Jona-Lasinio (NJL) model to describe its structure. The NJL model is not a toy model. In fact, it 
is the most realistic model for the pion based on a quantum field theory built with quarks. It gives a good description 
of the low energy physics of the pion and respects the realization of chiral symmetry . Moreover it has been used 
as the model to tune many coefficients of Ghiral Perturbation Theory [T^ . 

The NJL model is a non renormalizable field theory and therefore a cut-off procedure has to be defined. We have 
chosen the Pauli-Villars regularization procedure because it respects all the symmetries of the problem. In our scheme, 
we use the NJL model to describe the soft (non perturbative) part of the process, i.e. the initial and final states, while 
for the hard part we use conventional perturbative QCD. The use of the NJL model allows to calculate the GPDs for 
massive pions. Some peculiarities, as the non vanishing of the GPDs at the kinematic boundary regions, x = 0 and 
X = 1, well known in the exact TOtt = 0 limit, survive when m^r 0. 
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FIG. 1: Diagrams up to order g for the GPDs. The black dot indicates the scalar interaction vertex, while the white blob 
indicates the effective quark - gauge-boson vertex. 


Our paper is organized as follows. In section II we give some general definitions for scalar partons and we introduce 
our kinematical variables. In section III and IV we will define our approach for the scalar and the NJL model 
respectively. Section V presents our results and Section VI our conclusions. 


II. GENERALIZED PARTON DISTRIBUTIONS 


The GPDs are non-diagonal matrix elements of bi-local field operators. Various conventions, reference frames, vari¬ 
ables, etc., have been used in the literature for the description of such objects. Our notation is represented in Fig.^ 
i.e., the initial momentum is labeled by P, the final momentum by P', and the momentum transfer is given by 
A = P' — P. We shall describe initially a model with scalar particles, the generalization to particles with spin is 
straightforward and will be outlined in a later section. The GPD of a scalar system is defined by the matrix elements 
of bi-local scalar field operators Hiii: 


J+ = 


1 

2 


„ixP+z~ 

2tt 


(0)5+$(z)|P) 


z+=z-'-=0 




( 1 ) 


where x is the conventional Bjorken variable, ( the so-called skewedness parameter, and d = d — d. The elastic 
electromagnetic form factor of a system composed of two scalar particles is given by: 


j+ = (p'l $t ( 0 ) 9 + $ ( 0 ) \p) = {P + P')+ F{t). 


( 2 ) 


It follows directly from these definitions that integrating the GPD over x gives the form factor, 


J C, t) dx = F{t), (3) 

where the dependence on the skewedness parameter C drops out. This result is an important constraint for any model 
calculation. 

Let us now present our notation. Any four vector will be denoted (u"*",u"*",u“), where the light cone variables 
are defined by = (u° ± v^)/'/2 and the transverse part v-^ = For the kinematics indicated in Fig. ^ we 

introduce the ratios 


X = 


P+ 



( 4 ) 


of plus-components. With these definitions, which differ from other conventional ones [HQ , both X and C are 
defined on the interval [0,1]. 

We are only going to consider elastic processes, so P^ = P'^ = and = t. The following relation is true in 
general: 


(A-^+CP^)'=-C"M2-(l-C)t, (5) 

the positivity of which implies an upper bound for the skewedness ^ at a given value of the momentum transfer t: 

C < (-t)/(2M2) (v^l + 4M2/(-t)- l) < 1. 
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III. SCALAR ELECTRODYNAMICS 

We begin our study with the simplest model which allows for a completely analytic solution of the Bethe-Salpeter 
equation: this model describes a bound state of two distinguishable equal-mass scalar particles bound together by 
a zero-range interaction. For later convenience, we choose only one of the constituent particles to be charged. Our 
Lagrangian is, 


£ = [D^(j)]^ [D^^cj)] - - | , ( 6 ) 

with + ieAfj, so that the electromagnetic charge only couples to the field (j). Assuming that the coupling 

constant g is larger than a critical value, we have bound states stemming from the last term in Eq. The 

corresponding Bethe-Salpeter equation is trivially solved in the ladder approximation M , providing the fully covariant 
amplitude for the bound state of total mass = P^: 


d>(P,p) = 


C 


(jp — m? + ie) (^{P — — w? + 


(7) 


where p is the four-momentum of one of the two constituents with equal mass m, and C is a normalization constant. 
The spectrum for the ground state in this model, can be obtained directly from the Bethe-Salpeter equation, and is 
given by 


1 = -ighiP), ( 8 ) 

where the integral hiP) is given by Eq. llB^ . Note that this integral is actually divergent, so some sort of renormal¬ 
ization would be required in order to calculate the spectrum. The theory defined by Eq. is renormalizable and a 
renormalization program for bound states can be defined. However, for the evaluation of parton distribution functions 
in this model, we do not encounter any divergent integrals, we shall therefore ignore any matters of renormalization. 

The model defined by the above equations is certainly not realistic enough to furnish a reasonable description of real, 
physical bound states, like the pion for instance. Its main advantage lies in its simplicity, the fact that one may obtain 
analytic solutions, avoiding approximations that might destroy physical requirements or symmetries, like Lorentz- or 
Gauge invariance, sum rules, etc. Moreover it has the added advantage that it is defined within a renormalizable 
quantum field theory. These properties make it a useful playground to perform benchmark calculations, as it was 
used recently in order to test the viability of certain relativistic quantum mechanics approaches [laill. 


A. Generalized parton distribution function 

Our starting point is to write the generalized parton distribution (GPD) as an integral over Bethe-Salpeter amplitudes. 
To do so we could use Eq. o, a procedure which we shall develop in the NJL case. However for completeness we 
follow here the procedure of ref. 1^ . which establishes a relation between the calculation of the GPD and that of 
the electromagnetic form factor. Namely, in the scattering process 7 * -I-tt —> y-l-Tr, the leading “hand-bag” diagram, in 
the deep inelastic limit, reduces to a triangle diagram with an effective vertex (containing two-photon contributions), 
resulting from the contraction of the propagator with infinite momentum (see first diagram in Eig. ^|. Then, 

n{x,C,t) = J + [(-P-p)^-w^] ^’(-P + A,p-bA) 

Onp _ 1 r — 

= —^—- J - P'^/P^)<^{P,p)[iP - P)'^ - + ^,P +^)- ( 9 ) 

Note that the adjoint Bethe-Salpeter amplitude $ is defined via an anti-chronological time ordering of field opera¬ 
tors |l 9|| as compared to $, which means that the infinitesimal part ie in Eq. 0 keeps the same sign as compared to 
the mass term m^. 

Using the Bethe-Salpeter amplitude of Eq. 0 , we obtain 

2x — ( C'^ f d*p S{x — p'^/P~^) 

(27r)"^ (p^ — m'^ + ie) [(p -f A)^ — m'^ + ie] [(P — p)^ — -I- ie] 


n{xX,t) 


2 


I 


( 10 ) 
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Inspecting the pole structure of the integrand for the evaluation of the p integral, we note that it vanishes unless 
0 < a: < 1, i.e., the GPDs have the correct support properties. The integral of Eq. ns may be calculated explicitly. 
The analytic result for the GPD is then simply 

(j 2 2x — C~ 

0 < a: < 1, (11) 

with l3{m,xX,t) given by Eq. ifRbl . We observe that our result is explicitly covariant, depending only on x, C 
and t. Note that the complete solution is explicitly continuous at x = ^, however, the derivative at this point is 
discontinuous. We furthermore encounter a zero at the point x = C/2, which is due to the photon vertex, see Eq. I©. 
The quark distribution function is given by 


q{x) = Ti.{x, 0 , 0 ) 


x(l — x) 

IGtt^ — x(l — x)M^ 


( 12 ) 


The normalization integral may be done analytically and determines the normalization constant C. 

As it is not a common practice to write GPDs as integrals over Bethe-Salpeter amplitudes, we note that we can 
write 7f(x,C,t) for x > C as the product of light cone wave functions, defined by: 


4'(x,p-^) = ^ — f dp $(P,p) = -- 

ITT J I 


C^Jx{\ — x) 


(p-*- — xP-*-) + w? — x(l — x)M‘^ 
which is non-zero only for 0 < x < 1. Using some kinematic relations, the GPD for x > C may then be written as 


(13) 


'H{x,C,t) 


x>c 


1 2x- C 

2 2 




d‘^p 


\/x{x - C) J (27r)' 


-vl/* 


x-C 

1 -C’^ 


1 — X 


+ ^A^)vI;(x,p^). 


(14) 


On the other hand, for x < C, tbe GPD may not directly be written as the product of light front wave functions like 
in Eq. da, even though we are still able to find analytic solutions for the integral of Eq. m- 

Another conventional way of writing GPDs is via a parameterization in form of double distributions 0. The 
relation of these double distributions and the associated D-term with the results of a scalar model similar to the one 
considered here was recently discussed in ref. |^. 

In Fig.|2we give examples of GPDs in this model for different values of the binding energy and momentum transfers. 
We distinguish two scenarios: i) the weak binding scenario M « 2m (up to small binding energies of electromagnetic 
magnitude); ii) the deep binding scenario, which we choose close to the chiral limit, M = « 140 MeV, for reasons 

which will become natural in our later study of the Nambu-Jona-Lasinio model. 

We checked numerically that the sum rule of Eq. m is always exactly satisfied. A more stringent test is the so-called 
polynomiality condition [J, which states that the moments of the GPDs, 

[ dx x^~'^n{x, C, t) = Fn{C, t), (15) 

Jo 

give functions FAr(/',f) that are polynomials in ^ of order not higher than N. It is difficult to verify this analytically 
in the general case, even with the exact solutions that we obtained. Just in the limiting case = 0, the following 
relation may be shown to hold: 


Fw(C,0)=3 


2N+{N-2)j:tiC 

N{N + 1){N + 2) 


(16) 


i.e., the polynomiality condition is exactly satisfied at zero momentum transfer. Note also that for N = 1, Eq. itTCI) 
gives the correct normalization of Eq. © for any value of C. This shows explicitly that the sum rule is indeed 
independent of / in this case. For t ^ 0 and M'^ 0 we have shown, by numerical integration, that the polynomiality 

condition also holds (see appendix IHll. 

In the small binding limit, ^ 4m^, we recover the peaked nature of the GPDs, that was observed in ref. [^ . 
where also an interpretation was given. For small binding energies, the constituents are almost free, so their wave 
function is highly peaked in momentum space. The two peaks at x = 1/2 and x = (1 -I- C)/2 correspond just to the 
maxima of the corresponding wave function in the overlap formula Eq. OH- In the exact limit M = 2m, the GPDs 
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FIG. 2: Generalized parton distributions in the scalar model for different values of the bound state mass M and momentum 
transfers t. The top row gives results for t = —0.1 GeV^. The lower row those for t = —1 GeV^. The graphs follow the 
scenarios discussed in the text: i) (weak binding, M « 2m) left column and ii) (deep binding, M — m-„) right column. The 
values of the third variable are shown in the figures. Note that the sum rule of Eq. 0 is exactly satisfied for each graph. 


reduce to a sum of two d-functions^: 






(17) 


In particular, the quark distribution function q{x) = 'H{x, 0, 0) becomes simply q{x) = 5{x — 1/2), which means that 
the particles are free. 

In the deep binding limit, M = the GPDs change from the sharp peaks of the previous case to broad bumps 
describing the Fermi motion of the deeply bound constituents. At zero momentum transfer, t = 0, we find the simple 
form 


7f(a;,C,0) = 6^^^^ |^6»(C- a:) 6 [x) + ( 2 ; - C) (1 - f- • 


(18) 


B. Electromagnetic form factor 

Calculating the electromagnetic form factor according to Eq. we find: 

d'^p (A + 2pY 


[P + P'YF{Q^)=iC^ J 


(27r)4 (p2 _ ^2 _|_ _ pY _ ^2 _|_ jgj ’ 


(19) 


Note that the form factor in this limit is equal to zero everywhere except at t = 0, where it is 1. 


1 
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FIG. 3: The second diagram of Fig. ^contains the loop diagram shown on the left, which corresponds to the expanded set of 
diagrams on the right, including the seagull term. All of them are divergent, but their sum is not. Moreover when we sum them 
exactly also there hnite parts vanish. The reason behind this cancellation is that the charge is not renormalized to order g. 


leading to 


F{Q^) 


f\ 2 1-z, VD + zQ 
Jo Q Vd ^^Vd-zQ 


2-C, 


dxT-L{x, 


( 20 ) 


with D = A[m? — z(l — z)M^] + z^Q"^. This suggests that one could calculate the GPD by a simple change of the 
integration variable in the expression for the form factor. The general character of this change of variable would not 
be clear, however. Naturally, at = 0 we reproduce the normalization condition. 

Taking the parameters m = 241 MeV and M = 139 MeV, we can evaluate the root-mean squared radius of a “pion” 
built of scalar particles, via (r,r) = —6 = (0.47fm)^, to be compared to the experimental value 

of (0.66fm)^. In the limit of a mass-less bound state, M = 0, we obtain the analytic result (r,r)^ = 3/(10= 
(0.45 fm)^. 


C. Bubble diagram and Gauge invariance 

Up to first order in the strong coupling constant g, we have to calculate the second Feynman diagram of Fig.^ In 
scalar Electrodynamics, the second diagram does not contribute to the electromagnetic form factor. The loop integral 
is proportional to 




(2g-h A)^ 


(27r)'^ (g2 — -I- it) {{q + A)^ — + ie) ’ 


( 21 ) 


which by virtue of Eq. 1B4II is zero. However, this diagram could contribute to the GPDs, because here we keep the 
-|- component of the integration variable fixed. In a naive way we observe that the bubble present in the diagram 
will lead to divergent contributions to the GPDs. In order to understand properly the problem let us consider the 
DVCS on the pion. Unfolding the procedure introduced in this section, we note that the bubble corresponds to a 
triangle in which one of the propagators has been contracted due to the infinite momentum carried by the parton. 
Divergences appear in the calculation and correspond in practice to terms proportional to log((5^) in the triangle 
diagram. Considering all contributions to the same order (see Fig. O, we find that the divergences cancel. But not 
only the divergent terms cancel, also the finite parts do, therefore the whole contribution, contrary to expectations, 
vanishes exactly for the GPDs. 

It is not difficult to advance a physical argument explaining this cancellation. The diagrams of Fig. |21 can be 
interpreted as a renormalization of the bare charge present in the seagull term. Due to gauge symmetry this charge 
must be the square of the charge associated to the one photon vertex, before and after renormalization. But we have 
seen from Eq. that there is no renormalization (finite or infinite) at order g to the charge. Consequently the sum 
of diagrams of Fig. 13 or the second diagram of Fig.Q] must vanish. 


IV. THE MODEL OF NAMBU AND JONA-LASINIO 

The model of the last section can certainly not provide a very realistic description of electromagnetic properties of 
a pion because of the scalar nature of the constituents involved. We are therefore going to consider a model with 
spinor particles. An evident choice to investigate is the model of Nambu and Jona-Lasinio (NJL) [n m m m. 
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This model is the most realistic model for pions built of quarks that is based on a quantum field theory. It gives a 
good description of the low energy physics of the pion exhibiting the phenomena of dynamical mass generation and 
spontaneous breaking of chiral symmetry, which are crucial ingredients for low-energy hadronic physics HEES 
We start from the Lagrangian density in the two-flavor version of this model where we add the electromagnetic 
interaction in the usual way, 


C = -f g 





( 22 ) 


with = d/j, + ieA^. Again, assuming the contact interaction of the last term to be mainly responsible for the 
binding, the Bethe-Salpeter equation in ladder approximation to be fulfilled by a bound state in this model is given 
by (the factor 2 comes from the symmetry of the interaction): 

S-\p)^{P,p)S-\p-P) = igT^^ I 2Tr{f75-$(P,p')}^- (23) 

Here, the symbol Tr refers to traces on spinor, flavor and color indices, and S{p) is the single quark propagator 


Sip) = 


i 

p — m + ie 


i 


p + m 

p^ — m? it ’ 


(24) 


of a quark with constituent mass m, which is generated from the bare mass /io via the gap equation Eq. 

The quantity <i>(P, p) is the momentum space image of the Bethe-Salpeter amplitude of a bound state with total 
four-momentum P, while p is the four-momentum of one of the constituents. The solution of Eq. 03) is rather trivial 
since the integral on the right hand side is just a constant, so we can write 


$(P,p) = igjrqqS{p)Tj5S{p - P)= -igTzqqJ^ 


{p + m)fj 5 {p — ]p+ m) 
w? -\- *e)[(P — p)^ — rrP + ie] ’ 


(25) 


where gT^qq is the quark-pion coupling constant (given in Eq. (IBllIl l which can be determined from the usual Bethe- 
Salpeter normalization condition . Reinserting this solution into Eq. 03 gives a self-consistency condition, (given 
in Eq. (IHToli l. which determines the mass of the ground state as a function of the coupling constant. Note that in 
the chiral limit, when P^ = 0 and po = 0, Eq. fTTmli is nothing but the gap equation IBQII . providing some evidence 
for the self-consistency of the procedure. 

The NJL model has been investigated quite extensively in different domains of physics with rather impressive 
success (see refs. Eii for reviews), it is therefore natural to test its predictions for GPDs of the pion. Similar 
studies have been carried out recently 11 m El El 113 with different models. One drawback of the NJL model 
is of course its non-renormalizability, which makes it useful only as an effective, low-energy model that, however, 
may be regarded as a non-linear realization of the QCD Lagrangian. Numerical results therefore usually depend on 
the regularization scheme employed to deal with the divergent integrals. As thoroughly discussed in refs. a 

suitable regularization method has to satisfy a certain number of requirements. The method that was found to be 
most suitable was a Pauli-Villars with two subtractions, this is the one that we shall adopt, as outlined in appendix IXI 


A. Generalized parton distribution 

We are interested in calculating the diagrams of Fig. 0] in the Nambu-Jona-Lasinio Model. A general expression for 
the matrix element of the bi-local M-quark current in the Bethe-Salpeter approach is 


(P'|4'„(a;')7^4'„(a:)|P) = / d'^a;2Tr $p,(a;', a;2)i(l-k T3)7^ - m2) 




-TOi)$p(a:i,a;')^(l + T3)7^|, (26) 


where i is the isospin index. The first (second) term in this expression corresponds to the contribution of the first 
(second) constituent of the system, and indices 1 and 2 refer to coordinates or operators related to particles 1 and 2, 
bound in the meson. In terms of momentum variables we have 


(xi,X 2 ) = e 


-iPX 


(Pk 


e-*'=’'4>(fc,P), 


(27) 











p 


p + A 


P 


p — A 


P 


p — P 


P’ P 


P +p 


P’ 


FIG. 4: Diagrams contributing to the GPD in the Nambu-Jona-Lasinio Model. 


with the center-of-mass and relative coordinates defined by 

X = HlXl + ^ 22 ^ 2 , r = Xl— X2, 


Ml ,2 = 


Wl ,2 

mi + 7712 ’ 


(28) 


and P = Pi + P 2 , k = P 2 P 1 — M 1 P 2 the total and relative four-momentum of the system. 

We shall concentrate on the u-quark parton distribution in a 7r+ meson, for which we have to consider the first 
diagram of Fig. 01 with pi = p,p 2 = P—p. The second diagram of Fig.^contributes to the d-quark parton distribution 
with Pi = —p, p 2 = P + p. Inserting Eq. m and taking mi = m 2 = m, the diagrams of Fig. 3] give the following 
contribution to the GPDs: 


1 

2 


dz 




= T-i {x,C,t) = 


z+=z-'-=0 




(29) 


where we have switched back to particle momenta. With the Bethe-Salpeter amplitude of Eq. ll^ we obtain an 
expression for the GPD in the Nambu-Jona-Lasinio model that is similar to Eq. (US: 


7f(x,C,i) = i4:Ncglqq 


<Pp _ + ip+\ ( 2 ^ + 1 - C)(P^ -rrp) Pp- X-xP ■ X -{2x-O p-P 
{2Try ^ ^ {p'^ — m'^ + ie) [{p + Ap — m'^ + ie] [{P — pP — m'^ + ie] 


(30) 


The p integral in Eq. (EOll is evaluated by the usual residue calculus. Due to the pole structure of the integrand 
we obtain two contributions, the first one in the region C < 2 ; < 1, corresponding to the quark contribution and the 
second in the region 0 < 2 ; < ^, corresponding to a quark-anti-quark contribution. The second diagram of Fig. ^ 
will give an anti-quark contribution in the region C — 1 < 2 ; < 0, and a quark-antiquark contribution in the region 
0 < X < just like the first diagram. In case of a tt'*" meson, the first diagram gives the it-quark distribution, 
while the second one gives the d-quark distribution, but due to isospin symmetry, both distributions are related 
by Hu(x,C,t) = —7dj(C — x,C,,t). Concentrating only on the rt-quark distribution and employing a Pauli-Villars 
regularization, we get for C < a: < 1: 


'H(x,C,t) 


C<a:<l 



m“ 1 , m^ — xXP 1 , m^ — yM"^ 


- log —I - - log 
2 


m] 


-2 log- 


ml; 


{2x — C)M^ + (1 — x)t 


h(,rnj,x,C,t) 


(31) 


with the abbreviations x = x{l — x), y = {1 — x){x — C)/(l ~ C)^? while l 3 {m,x,(^,t) is given by Eq. frm with 
m = mo. 

Turning our attention to the non-valence region, 0 < x < ^, we find a first contribution arising from the diagram 
of Fig. 01 which is given by 


7f“(x,C,t) 


Ncg: 


2 

T:qq 


0<x<C 


47r^ 


E< 

J=0 


X nP I riP — 

- - log - log - 2 - 

4 m^ 2 m^ 


2x/C - 1 - yt 


log- 


-b 


(2x — ()M'^ + (1 ~ x)t 


h{mj,x,C,t) 


2 


(32) 
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FIG. 5: Diagrams contributing to the GPDs in the Nambu-Jona-Lasinio model coming from the ”a” coupling to the two 
photons. They contribute in the region 0 < a; < (". 


where now y = x{(^ — x)/C'^. In the region 0 < x < there is a second contribution coming from the re-scattering 
diagrams of Fig. O They correspond to the coupling of the two photons in a channel with the quantum numbers of 
the a meson. Although not apparent, their contribution is of the same order as the one from the diagrams in Fig. 0] 
(see appendix El for details of the calculation). The result for the sum of these two diagrams is 


H^{x,C,t) 


Ncg: 


2 

Tvqq 


0<x<C 


47r^ 


2 


^ 1=0 


- In - In „ — 

mi 


(33) 


with 

h{P,P') + 2h{P-P')] 
^ ^ (M2) -b (4 to2 - t) h (P - P') 


(34) 


and y = x{( — x)/(‘^. 

In the exact chiral limit, M = ttIt^ = 0 and rria = 2?7i, Eq. 121 becomes 


C{t) 




h{P,P') ' 
hiP-P')^ 


(35) 


where the tr propagator appears in an explicit way. This propagator is also present in Eq. but it is not so 

apparent due to the more complicated structure of the equation. Thus we can speak effectively of a coupling of the 
two photons to the a state, and this contribution is a requirement of chiral symmetry. 

The total value for the GPD in the region 0 < a; < C is obtained by summing the two, 


n{x,Q,t) 


0<a;<C 


n°-{x,c,,t) 


0<a;<C 


0<a;<C 


(36) 


We note that the first three terms in curly brackets of Eq. 121 give contributions that are independent of the 
momentum transfer t. The first term, in particular, gives an overall constant, independent of x, C, and t, which does 
not vanish for x = 1. From Eq. 121 we get the distribution function at zero momentum transfer. 


q{x) = H{x, 0,0) = 


Ncglqq 

47r2 



log ^ - log ■ 




m] 



(37) 


which is consistent with the normalization condition. 

In Fig. Ewe give some examples of GPDs in the NJL model for different values of the binding energy and momentum 
transfers. We distinguish in here three scenarios: i) the weak binding scenario M « 2m (up to small binding energies 
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of electromagnetic magnitude) ; ii) the deep binding scenario, with binding energies in the tens of MeV range, which 
would correspond physically to a strongly bound state, i.e. like the p meson; iii) the Goldstone boson scenario, which 
we take close to the chiral limit, M = « 140 MeV. This classification is very natural as the behavior of the GPDs 

in the figures shows. The Goldstone boson scenario, with its flat GPDs, did not appear in scalar Electrodynamics 
and it shows how spontaneous chiral symmetry breaking appears in the GPDs. 

For zero binding, there is no spontaneous symmetry breaking, thus no pion, and the constituent quarks are quasi 
free. For weak binding, we expect for the GPDs a similar behavior to the one discussed in scalar Electrodynamics. 
This similarity can be seen by comparing the corresponding graphs of Figs. [3 and In the weak binding limit we 
note again the appearance of two peaks at a; = 1/2 and a; = (1 + C)/2, which go over into J-functions in the exact 
limit M = 2m, just as in the scalar case, see Eq. EJ. As soon as there is binding the GPDs differ more and more 
from the free ones as we approach the chiral M = limit with massive constituent quarks. 

A few properties of the final result may be recognized from the analytic expressions of the GPDs. First, the total 
GPDs are discontinuous at x = (/. The discontinuity arises from the diagrams in Fig. [3 since the contribution from 
the diagrams in Fig. 0 to the GPDs is continuous, even though the derivative at x = C is not. Second, we see that 

(0, (,t) ^ 0 for any value of Its value, which depends on t and M but not on (, is connected to the discontinuity 
by: 

n (0, ct) = n (C+, C, t)-n (C- .C,t)y^o (38) 

Third, as noted above, H(l, t) ^ 0. Its value is independent of t and (, and only depends on the bound state mass 
M via the coupling constant gT^qq. In particular, in the chiral limit, when M = 0, the value of 7i(x,/',t) at x = I is 
simply given by: 

n{i,c,t) = i. (39) 

Furthermore, the distribution functions become relatively more concentrated around x = I for large |t|, which is just 
a consequence of the constant value of the distribution function at x = 1. Finally, in spite of these peculiarities, we 
find again that Eq. © is always exactly satisfied, as we checked numerically. 

The polynomiality condition in this case takes a very simple form in the chiral limit. When M = 0, the following 
relation holds: 

F^(C,0) = 1 (^I-ic"), (40) 

i.e., at zero momentum transfer, the polynomiality condition is fulfilled. Again, for N = 1, Eq. (gnj) reproduces Eq. © 
in the case t = 0. For t ^ 0 and ^ 0 we have shown, by numerical integration, that the polynomiality condition 
also holds. Further details are given in appendix I dI 

The fact that the polynomiality condition holds both for SED and NJL in the present scheme should not be a 
surprise since both models arise from a field theoretic description which preserves Lorentz symmetry, parity, time- 
reversal invariance and hermiticity 0. 

Finally we note that in the chiral limit, M = 0 and at zero momentum transfer, t = 0, we reproduce the analytic 
results of ref. 0 , namely: 

n{x, C, 0) = (C - x) 6» (x) -f 6» (x - C) 0 (1 - x), (41) 

so in particular, 7f(x, 1,0) = 1/2 and 7f(x,0,0) = 1. It may be checked from Eq. (T^ that this result is completely 
independent of any regularization procedure. This form may be compared with the corresponding result in the scalar 
model, see Eq. (UHl). Eq. (gU also implies that in the chiral limit we get a quark distribution function that is equal 
to unity. This coincides with the results obtained in refs. [^1^1^. 


B. Isospin decomposition 


Following the idea of ref. we introduce an isospin decomposition for the GPDs of the pion^ 


dz 




27T 


^ Note that our notation is slightly different from the cited reference. 
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FIG. 6: Generalized parton distributions in the Nambu-Jona-Lasinio model for different values of the bound state mass M 
and momentum transfers t. The top row gives results for t = —0.1 GeV^, while on the bottom for t = —1 GeV^. We describe 
from left to right the scenarios discussed in the text, i.e., from weak binding {M « 2m, left column), to strong binding (middle 
column) to the chiral limit {M = 0, right column). Note that the sum rule of Eq. 0 is exactly satisfied for each graph. 


In order to extract the two isospin components we need to calculate the rt-quark GPD for the 7r°. We have contributions 
coming from both diagrams of Fig. 0 and from those of Fig.0 leading to 


- niC - xX,t)], 


(43) 


where H (x, C, t) is the u-quark GPD for the 7r+ defined in Eqs. (13111 and II36II . Here x ranges from ^ — 1 to 1, as has 
been explained below Eq. (031. The GPDs of the 7r+ and 7r° allow to extract the isospin components 


°ix,c,t) =2nz ix,c,t) =nixX,t) -'Hic - xX,t) 

{xX,t) =2n^^ {x,c,t)-2nz° {xX,t) =n{x,c,t) + ?{{(- x,c,t) ■ 

Using Eq. 611 we obtain simple expressions for the two isospin GPDs in the t = 0, m.^. = 0 case: 


n 


1=0 


-1, C-l<a;<0 

{x,C,t) = { 0, 0 < a: < C 

1 , C < a; < 1 


{xX,t) = I, C~l<a:<l. 


(44) 


(45) 


We observe that the predictions of ref. [^, regarding the isospin components in the limit m,r = 0, are satisfied. The 
relation between {xX = Ij f = 0) and the wave function of the pion, (pTr{x), and (a;, C = 1, t = 0) = 0 hold. 
Note moreover that our pion wave function, in this limit, is ip^ix) = 1. The ttItt = 0, ^ = 1 and t = 0 case corresponds 

to a process in which the initial and final pion momenta are P = 0"*“: o) and P' = ^0, 0"'‘, 0^ , respectively. For 

zero pion mass this process is allowed only as a limit. In the NJL model we can look for deviations from the exact 
chiral limit. 


C. Electromagnetic form factor 

In order to understand better the above results, we may want to calculate the electromagnetic form factor in the 
NJL model which is given by the cc-integrated parton distribution, see Eq. This has been done already in several 
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works Ellllll, but we want to investigate the relation with the present results. Calculating again the current we 
find, from the first diagram of Fig.Q] for the form factor and for the GPDs: 


fP-L P'iM = S'At 2 [ + P^ + -m^+ ie)+p-APt^ - P ■ - (2p^ + A^)p ■ P 

[ + ) tv 1 o* "5-99 y (2^)4 (p2_^2+jg)[(p + A)2-m2 + ie][(P-p)2-TO2+ie] 

Rewriting the denominator under the integral and using some relations of appendix^ we arrive at 


(46) 


F{Q^) = t2N,g%^ 


{Q'^12 + M2)/2(A) + M^hiP) - M^hiA,-P) 


(47) 


We checked numerically (by applying the same Pauli-Villars regularization method) that this form factor is exactly 
reproduced by integrating our result for the GPD over x. This means that the non-vanishing distribution functions 
(at cc = 1) are implicitly present in the result of Eq. 14711 . Note that in the chiral limit, = 0, we simply have 

E(Q2) = 4tN,gl^ghiA) = 1 - RiQ^)/hiO), (48) 

where for small Q^, R{Q^) is given by 


R{Q^) 


167 r2 6rrP ’ 


(49) 


so together with the relation = — 12 im 2 j 2 ( 0 ), we can evaluate the root-mean squared radius of the pion in this 
model via (r^)^ = —6 9 F(( 5 ^)/ 9 Q 2 |^^^^ = which is the same as in ref. 0. With = 93 MeV, we find 

the numerical value of (r^)^ = (0.585fm)^, to be compared to the experimental value of (0.66fm)^. Note, however, 
that Eq. (1^ is valid only in the limit when the cut-off goes to infinity, for the finite cut-off as specified in appendix 1X1 
the root mean squared radius gets multiplied by a factor ~ 0.89. 


V. DISCUSSION 

The results obtained in the scalar Electrodynamics model show a perfect realization of all wishful ingredients. The 
calculation is exact, finite and satisfies all the desired properties, i.e., the GPDs have the correct support and vanish at 
the boundary regions, while the sum rule and the polynomiality condition are exactly verified. From the physical point 
of view, the GPDs show a realization in terms of quasi-free constituents in the weak binding limit. As the binding 
increases one is confronted with the dynamics as derivable from a non trivial momentum distribution determined by 
the corresponding Bethe-Salpeter amplitude, a feature also appearing in other model calculations 0 The physical 
effect associated with t and C is naturally represented in the GPDs. The variable t tends to push the constituent 
distribution towards higher values of a;, which corresponds to an input of momentum transfer into the system, while 
the variable ^ incorporates the description of virtual pairs. Unfortunately this model is not very realistic for the pion 
(perhaps it might be better fit for a description of the nucleon) and our results represent only qualitative features of 
how the dynamics might influence the distributions. 

An equivalent model was considered in refs. min, but the emphasis was put on the light-front quantization 
method. The results were calculated only for x > while the region x < (^ was explored by an analytic continuation 
of the vertex function. Due to the approximations used in their approach the continuation is not perfect and therefore 
sum rules like the one appearing in Eq. 0 are explicitly violated. 

For the NJL model the calculation requires regularization. The latter certainly influences the results, as has been 
discussed in ref. [sj, where the Pauli-Villars method was compared to the one of Brodsky-Lepage with different 
results. The Pauli-Villars method is compatible with all the symmetry requirements |33) . which is the reason for our 
choice. A caveat that should be emphasized is that the model does not only contain the dynamics expressed in the 
Lagrangian, but also the one derived from the regularization procedure. 

The NJL calculation retains some nice properties, in particular it preserves the sum rule and the polynomiality 
condition. Physically it is also very appealing since we can distinguish features associated with the weak and deep 
binding regimes. In Fig. [7|we show the variation of the quark distribution function with the binding energy. As 
the binding energy increases we change softly from a delta type behavior, at zero binding, to a constant behavior in 
the strong binding regime. This figure also illustrates nicely the phase transition associated with the spontaneous 
breaking of chiral symmetry. Suppose we keep the mass of the bound state fixed at = 0 and consider the variation 
of q{x) with the constituent mass m. It is clear that we will always have q{x) = 1 for any value of m, except when 
m = 0, where the distribution function changes discontinuously to q{x) = S{x — 1/2). The effects of t and ( described 





13 



- M^ = 0 


M^ = irrf 


- M^ = 3.9m‘ 


«- M^ = 3.99 m? 

|1 

1 1 

M 

1 1 

1 1 

I 1 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 


1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

[ ] 

-. 

^ . 


0 I"——. 

0 0.2 0.4 0.6 0.8 1 


X 


FIG. 7: Quark distribution function qix) = (*,0,0) in the Nambu-Jona-Lasinio model. The constituent mass is kept fixed at 
m = 241 MeV. 


for scalar Electrodynamics persist, as can be seen in the graph on the right hand side of this figure, i.e., t pushes the 
distribution to higher values of x and C introduces virtual pairs into the description of the system. 

Unexpected results of the calculation are the non-vanishing of the GPDs at the boundary, i.e., H{x,(,t) ^ 0 for 
X = 1 and a; = 0. This feature has nothing to do with what is usually called the support problem. The latter is 
characterized by a non-vanishing distribution function outside the physical region, i.e. ^ 0 for a; < ^ — 1 

or a; > 1, while in our calculations, the GPDs explicitly vanish there. It is therefore really a discontinuity that we 
encounter at the physical boundary. 

For the quark parton distribution q{x) = (x, 0, 0), this peculiarity was already noted in ref s. and it is 

obviously present in an implicit way in the results of electromagnetic form factor calculations [l^ ISft . It showed 
up already in several model calculations 0 , ■ 

Another peculiarity of the calculation is the discontinuity at the x = (. In the NJL model, the value of the GPD 
at X = 0 and the discontinuity at x = ^ are related, as we have shown in Eq. for any value of M, t and C- The 
value of H (0, C, t) decreases as t increases or as the binding energy decreases. But the origin of the non zero value of 
Ti ( 0 , (^,t) is in the cr diagrams depicted in Fig.|3 which originate from the scalar four-fermion coupling present in the 
Lagrangian, Eq. 

Let us discuss the discontinuity at the boundaries further, within the NJL description. One might investigate the 
mass dependence of the GPDs at a given point x. In our calculation we get for x = 1 the general expression 


n{x = i,c,t) 



(50) 


where Fq = 16171^/2(0), with l 2 {p) given by Eq. (IB6II and the dots denote small terms of higher order in 
where A is the cut-off parameter in the Pauli-Villars regularization, see appendix El This equation illustrates the 
dependence on the bound state mass, and also the regularization scheme dependence through Fq is apparent. It is 
clear from Eq. that 7/(x = 1, Cj^) vanishes for zero binding (M^ = 4m^) while it is non-vanishing, as soon as the 
interaction binds the quarks into a pion, i.e. < 4m^. 

Moreover in the NJL model, for a small mass particle like the pion (M = 

7/(x = 1, C, t) = + 0{mTr) -b ... = 1 -b 0{mTr) -b ..., (51) 


where the last step follows from the Goldberger-Treiman relation as a consequence of spontaneously broken chiral 
symmetry, and to lowest order in 77i,r this is regularization scheme independent. 

The physical interpretation of these features is not completely clear at present. Let us however conjecture a solution 
based on an analysis of known results. The usual idea that the parton distributions must vanish on the boundary 
is based in the study of free parton models. If one of the partons carries all the -b component of the momentum, 
i.e. X = 1, the other one has p'*' = 0 and thus p~ goes to infinity if the particle is on shell. Since this p~ enters 
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the propagator, the obvious conclusion is that the parton distribution must vanish at this boundary. But this simple 
description is broken when the particles are off-shell, like those in a bound system. Our result is therefore not unnatural 
and does not contradict any physical intuitive ideas. Moreover, the fact that the distributions do not vanish at a; = 0 
for C = 0 is just a consequence of the non-vanishing at cc = 1 and does not seem to introduce any new conceptual 
problem. 

One could suspect that this feature may originate as an artifact of the regularization in the NJL model, but this 
suspicion should be immediately dropped since it occurs in some limits in a regularization independent way. Therefore 
we conclude that there is a deeper physical reason for it, namely the off-shellness of the bound state quarks and the 
realization of chiral symmetry^. 

Up to now we have been discussing the NJL model as a theory in itself. However the true theory of the strong 
interactions is QCD and the NJL model should be interpreted as an effective theory of QCD at low energies. QCD 
based arguments led in refs. 0,113 to define a modified NJL model with a momentum dependent constituent quark 
mass (or equivalently, a constituent quark form factor) which, when appropriately chosen, makes the GPDs vanish at 
the end points. In ref. ^3 the emphasis was placed in the perturbative aspects, i.e., gluons and sea bremsstrahlung. 
The QCD evolution equations were used to eliminate the high momentum components. Certainly both mechanisms 
should be simultaneously advocated. 


VI. CONCLUSION 

In this work, we presented a detailed calculation of generalized parton distribution functions using a covariant Bethe 
- Salpeter approach both in scalar Electrodynamics and in the Nambu-Jona-Lasinio model. No assumptions have 
gone into the determination of the Bethe - Salpeter amplitudes or any other ingredients of the calculation. The only 
approximations employed are the ladder approximation for the determination of the Bethe - Salpeter bound state 
amplitude (which is, of course, still a fully covariant object in this case), and in the determination of the current matrix 
elements, we have restricted ourselves to the lowest order diagrams (but we have shown that the next order correction 
vanishes exactly). As a result of this procedure, no important features required by general physical considerations, 
like symmetry properties, sum rules, etc., have been violated and we recover them in our numerical results. For the 
scalar model, the calculation evidences all desirable features, we reproduce the results obtained in similar studies, 
extend them to other kinematical regions and find sum rules which were not possible in the other treatments. 

In the case of the NJL model, we found that the realization of chiral symmetry plays a crucial role in the outcome 
of the calculation. In the region 0 < a; < C chiral symmetry intervenes explicitly through the contribution from 
the scalar channel. In the region C < x < I, where the cr does not contribute, it appears implicitly through the 
consistency equations which govern the dynamics of the pion and its coupling to quarks. In the massive case, this 
reflects in the fact that the CPDs do not vanish at the boundary of the kinematic region, i.e., at a: = 0 and x = 1. 
Our detailed analysis associates this feature with the off-shellness of the constituent quarks, and the specific value 
there, with spontaneous chiral symmetry breaking. Moreover the CPDs are discontinuous at x = Q a consequence 
of the cr channel contribution,i.e., also a requirement of the chiral symmetry realization. These behaviors have been 
shown to arise in a regularization independent way. 

In conclusion, we have shown that the SED model has continuous, well behaved CPDs that vanish at the end-points 
of the kinematic region, and bubble diagrams which do not contribute as a consequence of gauge symmetry. On the 
other hand, the NJL model has discontinuous, ill behaved GPDs, that are non-zero at the end-points, and bubble 
diagrams which, in the a channel, do contribute as a consequence of the realization of chiral symmetry. 

Before finishing we must recall that our calculation is valid at the hadronic scale, i.e., at a low momentum renor¬ 
malization point. Evolution to higher momenta is necessary to describe deep inelastic scattering data. It will be 
interesting to see how the described features of the NJL model will change under evolution. However the fact that 
the distributions do not vanish at the boundary imply the appearance of strong singularities which render the process 
non trivial. 


^ Similar features occur in the well defined ’t Hooft model a theoretical scheme which is exact in the large Nc limit. In the chiral 
limit, this model produces a parton distribution q{x) which is equal to 1 for any value of ic, in a regularization independent way. 
But, when chiral symmetry is explicitly broken, the model of’t Hooft leads to a parton distribution which vanishes at the end-points, 
q(x = 0) = q(x = 1) = 0. This different behavior is not surprising since the dynamics of this model and that of the NJL model are very 
different. The’t Hooft model is a 1-hl dimensional model, which is confining, but without a point-like interaction. 
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APPENDIX A: REGULARIZATION 


We have used the Pauli-Villars regularization in the Nambu-Jona-Lasinio model in order to render the occurring 
integrals finite. This means that for integrals like the ones defined by Eqs. JbTIbI . we make the replacement 


/ 


d'^p 





(Al) 


with = mP + j'A^, cq = —Cil2 = C 2 = 1. Following ref. M we determine the regularization parameters A and m 
by calculating the pion decay constant and the quark condensate in the chiral limit via 


f 

J TT 


3m^ 

47r^ 


2 

^Cjlog(m2/m2), 

3=0 


(uu) 


3m 


2 

Cjin^ log(m|/m^). 

3=0 


(A2) 


With the conventional values (uu) = — (250MeV)^ and /tt = 93 MeV, we get m = 241 MeV and A = 859 MeV. 


APPENDIX B: ELEMENTARY INTEGRALS 


Some of the interesting integrals appearing in the main text are: 


f dH 1 
J {2'kY — m? +ie' 


(Bl) 


h{p) 


d'^k _ 1 _ 

(27r)4 (fc 2 _ 77^2 _|_ _|_ p'j2 _ 77 j 2 _|_ ’ 


(B2) 


h{pi,P2) = J 


d-^k 


(27r)^ (fc2 _ 77 j 2 _|_ jg^ _ 77 j2 _|_ ic)[[k +^ 2 )^ — m^ + ie) ’ 


From these definitions the following relation may be deduced: 

f d'^k k^ 


{2ttY {k"^ — rn?-\-ie) {{k pY — w?ie) 2^ 


(B3) 


(B4) 


The scalar model is renormalizable. The NJL model is not. We use in the latter case the Pauli-Villars regularization 
procedure and obtain: 


I^=- 


IGtt^ 


y]cj m| log(m|/m2). 


(B5) 


3=0 


hip) = 


i 

IGtt^ 


2 




log(m^/m^) 



(B6) 
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From Eq. IIB3II and p\ = P 2 = we have: 


hiPi,P2) = hit = (pi - P2f) = X! 

j=o 



^3 + 

5 

0 

Cj — \Jdj 


with 


dj = (to| + M'^z (z — 1)) , Cj = 2 to| + 2 M‘^z (z — 1 ) — z'^t. 

The self consistency relation determining the constituent mass of the quark, m, is 

m = po + 8 i g Nc Nf m Ii , 

while the equation for the mass of the pion , leads to 

l = 4.igN,Nf [2h - P^h (-P)] , 

and the pion quark coupling constant , g-Kqq, is given by 

2 _ 

9-Kqq 7 

12z (I 2 {Mh + M2 {dh ip) ldph^.=Mg 

To obtain the results of Eqs. inidllTlIl . we used the following formulas: 


hiP) = -16z7r2 


d'^k 


5ix-k+/P+) 


(27r)'^ (fc 2 — m? + ie) {{k — P)^ — m? + ie) I i=o 


I m'^^-xP^ 


0 


, 0 < a: < 1 , 

otherwise. 


(B7) 

(B 8 ) 

(B9) 

(BIO) 

(Bll) 


(B12) 


hixn,x,C,,t) = —16z7r^ 




Six - fc+/P+) 


(27r)"‘ (fc2 — vP + ie) ((fc + A)2 — m? + ie) ((P — k)^ — m? + ie) 


1 [ 2(1 - OiCy -x)- C]M^ - (1 - x)t + 2nP/y + /P q ^ 3 , ^ 

= <1 VP [2(1 - C)(C2/- a:) - C]M2 - (1 - x)t+ 2m2/2/-v/p’ ’ (3^3^ 

0 otherwise. 

In the last integral, D = C^M^(M^ — 4m^) + (1 — x)'^P + 2(1 — a;)(2a; — QM'^t — Aw? il — C)t, C = /P'^i M^ = P^, 

t = and 


_ x/(, 0 < X < h 

y ~ (1 - a;)/(l - C), C < a; < 1 . 


(B14) 


1 f d'^p ( p 


APPENDIX C: CONTRIBUTION OF THE a - DIAGRAMS 

The contribution of the diagrams shown in Fig. [3 can be calculated as T? (a;, (^,t) = A ■ B, with 

A = 

and 
B = 


I .Six— 4— 1 Tr 

2 7 (27r)^ V P+ 


iS ip) 7 +t (1 + r^) iS ip + A) 


‘^alqng 


d‘^k 


(Cl) 


i-25n,(p-P0 7 (2^)" 


Tr [ 75 T S ik — P) 'y 5 T~^S ik + A) S (fc) + (fc + P') 'y^r S ik + A) S (fc)] , 

(C2) 



















17 


where Us is the scalar vacuum polarization: 


ct^p 


n« (P) = -i / -^Tr [iS (p) iS {p - P)] = AiN.Nf 


h + - (4m" - P") h (P) 


After some algebra we obtain using Eq. (IHTHIl 


P}" (x, C, t) = iC-x)0 (x) (^1 - y ^ C (t) ^ Cj 


mi mi — ty 
_ In ^ - In ^ „ 
m^ m* 


with y = x(C — x)/C" and 


C (t) = m 


2 [{t - 2M2) /3 (P, P') + 2/2 (P - P')] 
M2 /2 (M2) + (4m2 - t) I 2 (P - P') 


which are Eqs. 1221) and 122J. 


(C3) 


(C4) 


(C5) 


APPENDIX D: POLYNOMIALITY. 


One stringent test of models is the polynomiality condition. This study appears simpler and more physical in terms 
of the more symmetric variables 


^ a^- C/2 
l-C/2 


C 


C 

2-C' 


(Dl) 


With this definition X ranges from —1 to +1, and not from C — 1 to +1, which is the range associated with x. In 
terms of these new variables the polynomiality condition is stated as 





2 = 0 


it) e 


(D2) 


where [...] means the integer part. Time reversal invariance and hermiticity imply that the distributions (CiO 
even functions of C 0|- 

For SED in the case m^ = 0 and t = 0, as function of the new variables, we have 


7/(A,C,0) 


Using Eq. (ina we obtain for the coefficients, 


6 ^(X+C) 

1 +? 2 C 
6 X(l-X) 
1 +? 1 -? 


-c < X < c 

C < a: < 1 


A 


SED 

71,22 


( 0 ) 


6 

(n + 1) (n + 2) ’ 


(D3) 


(D4) 


This simple closed form only appears for vanishing t and m,r- For non zero values we have to proceed numerically. In 
tabled we show values of (t) for the lowest values of n, and non vanishing t and m.^-. These results show that 

F) anymore independent of i. We also obtain that the effect of a physical pion mass is small, i.e. at the 

level of a few percent. 

The definition of the coefficients in Eq. iTDll leads to the following relation for their sum, which is connected with 
the value of the moments at C = 1, 


i/2m-l(l,0) =P2m(l,0), (D5) 

with m = 1,2,... We have verified numerically (with a precision of 10“^^) that this relation is also exactly satisfied 
for t 0 and = 0: 


771—1 


H2m-1 (1, t)=Y, it) = H2m (1, i) = ^ F) > 


2 = 0 


2=0 


(D6) 
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For 7 ^ 0 the maximum value of ^ never reaches 1 and therefore we cannot connect the sum of the An^ 2 i coefficients 
with the value of the moment at any physical value of Our numerical study shows however that, in the case of 
SED, we have an approximate sum rule for the sum of coefficients An^ 2 i' 


m — 1 

E 

2=0 


(t) - E {*) - ^ (to^) 


2=0 


(D7) 


The deviation from the exact sum rule is at the level of one percent for t = —10“^ GeV^ {—t small) and at the level 
of 10 “^ for t = —10 GeV^ (—t large). 

We now turn to the NJL model. As in the SED case, we start from the parton distribution for t = 0 and m,r = 0: 

n{x,^,o) = - X) e {^ + X) + e {X - 0 e{1 - X) (ds) 

Inserting Eq 1D8II in Eq. (ID2II and integrating over X we obtain the coefficients. 


< 0 ^ ( 0 ) 
AZiHo) 
( 0 ) 


1 

n 

0 , l< 2 z<n-l. 



(D9) 


In this case, the closed forms appear only for vanishing values of rriTr and t. Results for t 7 ^ 0 and TOtt 7 ^ 0 have 
been obtained numerically and are shown in table ||| From these results we realize that A^ 2 ^ (^) not zero for 
1 < 2i < n — 1. We also observe that the effect of a physical pion mass is small. 

We have found also several simple relations like 


i?2m(l,0) = 0, i72m+l(l,0)=i72m+l(0,0), (DIO) 

with m = 1, 2.... But in this case these relations are only approximately satisfied for t 7 ^ 0 and = 0 or for TOtt 7 ^ 0. 

We have studied numerically the contributions to (^, t) coming from each diagram separately. We observe that 
there are independent sum rules for each diagram but they disappear when the sum is performed. The contribution 
coming from the diagram in Fig. ^verifies the relation 

771 m 

E it) - E ^^+ 1 . 2 * it) = O {ml) (Dll) 

2=0 2=0 

exactly for m = 0 and any t, and in an approximate way (deviations less than 10“^) for m,n. = 140MeV. On the other 
hand, the diagram in Fig. Ogives non-vanishing contributions only to the coefficients A^l^ for any value of ttItt and 
any value of t. 

Regarding the numerical values of these coefficients we observe that in the scalar model and for large n, (0, 0) 
decreases with n as 1 /n^ whereas Hn ( 1 , 0 ) decreases with n as 1 /n, showing that Ti {X, 1 , 0 ) is more concentrated 
around X = 1 than Ti, {X, 0, 0). In the NJL case for large n, Hn (0, 0) and (1,0) decrease as 1/n. This shows that 
H {X, 0) is more concentrated around the X = 1 region than in the scalar case. These asymptotic behaviors appear 

for any value of t and m^^.. 

The numerical analysis also shows that for small n, as —t increases, the coefficients of the SED and NJL models 
are of the same order, loosing the characteristic behavior of equations and 1D9II . Nevertheless we observe a 
regularity in the signs: all the coefficients are positive in the scalar model, whereas in the NJL model the coefficients 
A„^q are positive but the rest are negative, except some small coefficients present at small —t. 

Finally we recall that Ai_o (0) = 1 corresponds to the charge sum rule and ^ 2,0 (0) = 1/2 is related to the momentum 
sum rule. 
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